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RAIN  CLUTTER  STATISTICS 


INTRODUCTION 

The  rain  clutter  data  presented  in  this  report  were  taken  with  a  frequency-agile  S-band  radar 
whose  parameters  are  shown  in  Table  1.  Five  equispaced  frequencies  from  2.1  to  3.9  GHz  were 
transmitted  and  were  varied  from  pulse  to  pulse  such  that  every  fifth  pulse  had  the  same  frequency. 
With  this  radar  a  single  antenna  is  used  for  both  transmission  and  reception  over  the  entire  range  of 
frequencies.  The  data  were  taken  on  5  September  1979  at  the  Randle  Cliff  facility  of  the  Naval 
Research  Laboratory  during  the  remnants  of  hurricane  David.  The  wind  at  the  surface  was  from  the 
northeast  at  approximately  23  knots.  The  facility  is  on  the  western  shore  of  the  Chesapeake  Bay,  and 
the  radar  coverage  lay  in  the  northeast  quadrant,  with  the  radar  looking  over  the  Bay  and  into  the  wind. 


Table  1  —  Radar  Parameters 


Frequencies: 

2.1,  2.55,  3.0,  3.45,  and  3.9  GHz,  varied  from  pulse  to  pulse 

PRF: 

1600/s 

Beam  width: 

1°  in  azimuth  and  2°  in  elevation  (at  3.0  GHz) 

Pulse  width: 

compressed  to  0.7  fis 

Antenna  rotation: 

0,  2,  and  6  RPM 

Polarization: 

vertical 

Our  main  concerns  were  the  functional  relationship  between  the  mean  backscattered  power  and 
frequency,  the  spatial  and  temporal  variation  of  the  backscattered  signal,  and  the  Doppler  power  spec¬ 
tra.  Among  our  results,  those  concerned  with  the  clutter  Doppler  spectra  have  the  greatest  potential 
importance  for  the  radar  engineer,  since  it  is  the  properties  of  these  spectra  which  determine  the 
effectiveness  of  clutter-rejection  filters.  A  Gaussian  model  for  rain  clutter  spectra  is  generally  accepted 
by  the  radar  community  (as  in  Ref.  1,  for  example),  and  the  early  experimental  works  of  Barlow  121 
and  Kerr  (3)  are  still  cited  as  primary  sources  on  this  subject.  However,  we  found  that  the  Doppler 
spectra  sometimes  had  multiple  peaks,  and  we  would  like  to  learn  something  about  the  nature  of  the 
processes  which  can  produce  such  multipeaked  spectra  and  how  frequently  they  can  occur.  The  sim¬ 
plest  explanation  for  the  occurrence  of  such  spectra  is  that  of  a  wind  shear  operating  on  a  vertically 
stratified  distribution  of  clutter,  but  there  are  other  possibilities.  One  such  possibility  frequently  dis¬ 
cussed  in  the  literature  is  that  of  Bragg  scattering  from  turbulent  eddies,  but  the  variation  of  backscat¬ 
tered  power  with  space,  time,  and  frequency  that  we  observed  is  more  consistent  with  the  vertical- 
stratification  model. 

From  the  spatial  and  temporal  variation  of  certain  statistics  we  found  that  the  observed  relation¬ 
ship  between  the  mean  backscattered  power  and  frequency  conformed  to  the  model  P(f)  «  AQif), 
where  Q  —  Q(f)  is  some  function  which  depends  only  on  frequency  and  A  is  a  frequency-independent 
parameter  which  varies  with  time  and  position.  For  the  Bragg-type  scattering  effects  discussed  in  the 
literature  we  have  Q(f)  —  f,  where  the  exponent  n  has  a  value  between  0  and  2.  But  from  the  data 
we  conclude  that  the  exponential  rate  of  growth  of  power  with  frequency  lies  in  the  range  3  <  n  <  7, 
with  the  preferred  value  being  approximately  n  —  5,  which  happens  to  be  the  value  of  n  which 
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corresponds  to  incoherent  scattering  from  narrow  horizontal  bands.  This  analysis  is  complicated  by  the 
possible  existence  of  calibration  errors  whose  effects  cannot  be  averaged  out,  because  they  are  constant 
and  multiplicative.  For  this  reason  (as  well  as  others)  we  did  not  estimate  the  value  of  n  by  fitting  a 
curve  to  the  data  in  the  usual  sense  of  least-squares  regression.  Instead  we  calculated  the  smallest  RMS 
calibration  error  that  would  be  required  to  produce  the  departure  of  measured  values  of  n  from  an 
assumed  "true"  value.  We  judge  that  the  RMS  calibration  error  was  no  larger  than  2  dB  (except  for  a 
glitch  at  the  highest  frequency,  to  be  described  in  the  section  on  system  calibration);  hence  values  of  n 
lying  in  the  range  0  <  n  <  2  are  rejected  as  reasonable  possibilities,  because  they  would  imply  an  RMS 
calibration  error  greater  than  3  dB. 

A  final  remark  concerning  methodology:  The  performance  of  radars  is  determined  by  the  state  of 
the  clutter  at  particular  times  and  particular  places.  Thus  our  interest  is  in  the  behavior  of  the  rain 
backscatter  over  short  intervals  of  time  and  localized  regions  of  space.  We  are  not  very  interested  in 
establishing  trends  in  data  which  are  heavily  averaged  over  space  and  time,  and  stationarity  was 
assumed  only  for  data  originating  from  a  single  range  cell  over  a  short  duration  of  time. 


ORGANIZATION  OF  THE  DATA 


The  basic  data  from  the  radar  consists  of  A/D  converted  /-  and  0-channel  voltages,  which  are 
organized  into  data  records,  each  consisting  of  28,130  16-bit  words.  These  data  records  are  temporarily 
stored  in  a  minicomputer  before  being  dumped  onto  tape.  Each  radar  pulse  gives  rise  to  a  block  of  data 
consisting  of  two  16-bit  words  recording  the  carrier  frequency  and  azimuth  and  Mc  16-bit  words  record¬ 
ing  the  /-  and  0-channel  voltages,  with  eight  bits  being  assigned  to  each  channel,  where  Mc  is  a  vari¬ 
able  indicating  the  total  number  of  range  cells  per  pulse  being  processed.  Hence,  the  total  number  of 
pulses  Np  in  a  record  is 


N, 


28,130 

2  +  M,' 


The  data  are  sampled  at  a  rate  of  3  MHz,  which  is  twice  the  reciprocal  of  the  compressed  pulse 
width  of  0.7  /is.  Hence  there  is  a  50%  overlap  between  successive  range  cells. 

When  the  antenna  is  rotating,  the  recording  device  is  triggered  to  take  data  centered  at  a  given 
azimuth.  Thus,  for  example,  there  is  a  10-s  delay  between  records  when  the  antenna  is  rotating  at  6 
RPM.  When  the  antenna  is  not  rotating,  there  is  approximately  a  2-s  delay  between  records. 

The  data  records  are  either  single-frequency  or  multifrequency  records.  The  single-frequency 
records  always  use  the  carrier  frequency  of  3.0  GHz,  whereas  the  multifrequency  records  use  each  of 
the  five  frequencies  given  in  Table  1.  In  the  multifrequency  case  the  frequencies  are  varied  in  order 
from  pulse  to  pulse;  hence  the  PRF  for  pulses  at  each  given  frequency  is  reduced  from  1600  to  320  in 
the  multifrequency  case. 

Finally,  for  convenience  the  data  records  are  grouped  into  files,  with  each  file  consisting  of 
records  taken  at  a  given  azimuth  and  centered  at  a  given  range.  In  effect  a  file  consists  of  data  from  a 
selected  clutter  cell  appearing  on  the  PPI  scope.  Table  2  gives  the  file  descriptions  for  the  files  referred 
to  in  this  report.  Of  special  interest  is  file  8,  which  was  obtained  by  turning  off  the  transmitter.  This 
file  will  be  referred  to  as  the  noise  file  and  is  used  to  measure  the  receiver  noise  statistics.  File  3  is  the 
only  single-frequency  file  used  in  this  report.  Figure  1  shows  a  schematic  PPI  display  of  the  regions 
covered  by  the  four  multifrequency  files:  1,  2,  5  and  6. 
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Table  2  —  File  Descriptions 


File 

No. 

General 

Description 

Antenna 
(RPM)  Rate 

Range 

(km) 

Az. 

(deg) 

El. 

(deg) 

No.  of  1 
Records  | 

Time  Duration 
per  Record  (s) 

No.  of 
Range  Cells 

~T~ 

Multifrequency 

6 

18.95 

359 

3.6 

0.27 

64 

Multifrequency 

6 

57 

3.6 

0.27 

64 

Single  Frequency 

6 

B&a 

57 

3.6 

5 

0.27 

64 

Multifrequency 

2 

31.78 

354 

3.2 

5 

1.76 

8 

Multifrequency 

0 

33.34 

40 

3.2 

1.76 

8 

LU 

Multifrequency  noise 

— 

— 

— 

— 

5 

0.52 

32 

NORTH 


VOLTAGE  AND  POWER  HISTOGRAMS 

Let  x„  and  y„  denote  the  /-  and  0-channel  voltages  measured  at  a  given  range  from  the  nth  pulse 
at  a  given  frequency.  From  general  principles  one  would  expect  x„  and  yn  to  behave  like  uncorrelated 
zero-mean  Gaussian  variables  with  the  same  variance: 

E[x„]  =  E\y„\  =  0, 

£[xn2)  =  E (>»„2]  -  (t2, 

E\xmyn )  =  0. 

Figure  2  shows  histograms  of  the  /-  and  0-channel  voltages  as  they  actually  come  off  the  A/D  con¬ 
verter.  Due  to  a  fault  in  the  A/D  conversion,  there  are  no  -  Is  and  Os  in  the  /  channel,  and  the  Os  and 
Is  are  underrepresented  in  the  0  channel.  The  main  effect  of  those  malfunctions  in  the  low-order  bits 
appears  to  be  the  production  of  small  systematic  biases.  These  biases  do  not  seem  to  vary  much  from 
file  to  file  (even  including  the  noise  file),  but  they  do  vary  slightly  with  frequency.  These  biases  are 
removed  numerically,  and  the  data  thus  obtained  are  said  to  be  bias  corrected. 

Table  3  shows  values  of  the  crosscorrelation  between  the  bias-corrected  /-  and  0-channel  voltages. 
The  values  shown  are  averages  obtained  from  all  range  cells  in  a  given  record,  at  each  of  the  two  indi¬ 
cated  files.  These  correlation  coefficients  are  small  in  magnitude  and  negative  in  sign,  but  they  are 
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VOLTAGE 

(a)  /-channel  voltage 


-12  -10  -8  -6  -4  -2  0  2  4  6  8  10  12  14 


VOLTAGE 

(b)  ^-channel  voltage 

Fig.  2  —  Voltage  histograms:  data  from  file  6  at  3.9  GHz 


Table  3  —  Crosscorrelation  Between  I-  and  Q-Channel  Voltages 


File 

Crosscorrelation  Coefficient 

1  2.1  GHz 

|  2.55  GHz 

1  3.0  GHz  1 

3.45  GHz 

1  3.9  GHz 
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larger  in  magnitude  when  signal  is  present  (the  coefficients  being  larger  for  file  6  than  for  file  8,  the 
noise  file). 

The  bias-corrected  data  were  used  to  generate  histograms  of  instantaneous  power,  and  these  are 
shown  in  Fig.  3.  Also  shown,  above  each  histogram,  is  a  plot  of  the  logarithms  (base  e)  of  the  number 
of  samples  in  each  class  interval.  These  plots  provide  a  check  of  the  validity  of  the  Gaussian  model  as 
follows.  The  data  are  from  multifrequency  files,  and  since  only  every  fifth  pulse  has  a  given  frequency, 
the  correlation  between  successive  pulses  at  a  given  frequency  is  small.  Hence  the  values  of  instantane¬ 
ous  power  for  each  pulse  at  a  given  frequency  should  be  independent  samples  of  a  variable  tv  =  x2  + 
y,  where  x  and  y  are  two  independent  zero-mean  Gaussian  variables  with  variances  both  equal  to  a2. 
The  variable  tv  has  the  probability  density  q  given  by 

q(w)  =  (l/w0)e  */w°,  (1) 

where 

tv0  =  E[w\  ~  2a2. 

The  expected  number  v  of  samples  which  fall  in  a  class  interval  of  length  A  centered  at  tv  is  approxi¬ 
mately 

v  =  NA<f(*v), 

where  N  is  the  total  number  of  samples.  Hence  In  v  should  be  approximately  linear  in  tv: 

In  v  —  -  (tv/tv0)  +  In  (AA/tVo). 

The  straight  lines  shown  in  these  plots  are  generated  by  this  relation,  with  w0  being  taken  to  be  the 
sample  mean  of  the  power.  As  can  be  seen,  the  fit  is  good. 

SYSTEM  CALIBRATION 

For  each  of  the  five  frequencies  used,  measurements  were  made  of  the  antenna  gain,  receiver 
gain,  transmitted  power,  and  cable  losses.  These  measurements  were  used  to  obtain  "calibration  fac¬ 
tors"  k,  (i  =  1,  2,  3,  4,  5)  which  are  proportional  to  the  reciprocals  of  the  total  system  response  at 
each  frequency.  Calibrated  values  of  power  were  obtained  by  multiplying  the  raw  values  by  the  factors 
k,.  These  "calibrated"  values  are  given  in  terms  of  unspecified  units  of  power,  since  no  attempt  was 
made  to  measure  the  power  in  watts.  Unfortunately  we  were  not  able  to  use  a  test  sphere  in  the  cali¬ 
bration. 

Previous  measurements  had  shown  that  the  transmitter  signal-to-noise  power  ratios  were  38  dB  or 
better  over  the  entire  frequency  range.  Hence  the  transmitter  noise  was  not  an  important  factor  in  our 
data  analysis,  since  the  corresponding  receiver  noise  levels  were  orders  of  magnitude  larger. 

Examples  of  the  calibrated  power  as  a  function  of  frequency  are  shown  in  Table  4  and  Fig.  4.  The 
table  shows  the  sample  mean  and  standard  deviation  for  each  of  the  indicated  range  cells  from  the  same 
indicated  record.  Also  shown  is  the  number  of  pulses  per  frequency.  The  power  is  calculated  from  the 
bias-corrected  data,  and  the  corresponding  means  and  standard  deviations  are  essentially  equal,  as 
would  be  expected  if  the  observed  values  of  instantaneous  power  were  random  samples  from  the 
exponential  distribution  (Eq.  (1)). 

The  power- versus-frequency  curve  shown  in  Fig.  4  was  obtained  by  averaging  all  the  data  in  file  6; 
however,  the  curves  obtained  from  data  in  any  individual  record  from  any  individual  range  cell  all  have 
the  same  general  appearance.  In  particular,  there  is  always  a  sharp  drop  off  in  power  at  the  highest  fre¬ 
quency,  3.9  GHz.  A  spectrum  analyzer  was  used  in  an  attempt  to  discover  malfunctions  in  the  signal 
synthesizer  or  transmitter  at  this  frequency,  but  none  was  detected.  Therefore  we  provisionally  ascribe 
this  anomally  to  an  unidentified  calibration  error.  As  will  be  seen,  the  existence  of  this  anomaly  has  no 
effect  on  our  conclusions;  in  particular,  the  returns  from  the  first  four  frequencies  were  sufficient  to 
exclude  Bragg  scattering  as  a  significant  factor  in  the  production  of  the  Doppler  spectral  peaks. 


5 


OCCURRENCES  NUMBER  OF  OCCURRENCES  NUMBER  OF  OCCURRENCES 


TTTiTi 


NRL  REPORT  8639 


Table  4  —  Mean  and  Standard-Deviation  Values  of  Power 
for  the  Record- 1  Data  From  File  6. 


Frequency 

(GHz) 

■UgjJ 

Bo 

Mean 

Std. 

Dev. 

Frequency 

(GHz) 

Mean 

Std. 

Dev. 

Range  Cell  1 

Range  Ceil  5 

2.1 

561 

35.2 

34.6 

2.1 

561 

28.7 

27.9 

2.55 

561 

104.4 

2.55 

561 

8,'  7 

80.8 

3.0 

561 

214.7 

211.4 

3.0 

561 

170.6 

176.6 

3.45 

561 

317.0 

323.7 

3.45 

561 

253.4 

259.1 

3.9 

562 

186.3 

182.4 

3.9 

562 

172.6 

171.9 

Range  Cell  2 

Range  Cell  6 

561 

31.8 

33.3 

2.1 

561 

32.4 

31.1 

■gN 

561 

97.1 

94.1 

2.55 

561 

89.3 

88.2 

3.0 

561 

214.1 

215.0 

3.0 

561 

159.5 

154.1 

3.45 

561 

284.8 

282.8 

3.45 

561 

255.1 

254.3 

3.9 

562 

167.2 

160.3 

3.9 

562 

180.2 

182.9 

Range  Cell  3 

Range  Cell  7 

2.1 

561 

29.3 

28.6 

2.1 

561 

33.0 

30.7 

2.55 

561 

96.2 

95.2 

2.55 

561 

88.5 

93.4 

3.0 

561 

210.5 

203.8 

3.0 

561 

179.6 

188.1 

3.45 

561 

274.5 

268.7 

3.45 

561 

264.6 

249.3 

3.9 

562 

176.8 

164.1 

3.9 

562 

178.0 

170.4 

Range  Cell  4 

Range  Cell  8 

2.1 

561 

28.2 

27.0 

2.1 

561 

30.7 

28.3 

2.55 

561 

91.0 

85.8 

2.55 

561 

91.6 

90.7 

3.0 

561 

194.8 

206.8 

3.0 

561 

199.3 

202.7 

3.45 

561 

263.1 

253.7 

3.45 

561 

309.1 

275.7 

3.9 

562 

176.9 

155.9 

3.9 

562 

1S8  0 

184.9 
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DOPPLER  POWER  SPECTRA 


Figures  5  and  6  show  Doppler  spectra  from  the  first  five  consecutive  records  from  the  third  range 
cell  in  multifrequency  file  6,  which  is  the  only  file  with  the  antenna  not  rotating.  As  stated  earlier, 
there  is  a  2-s  delay  between  consecutive  records  in  file  6.  The  PRF  for  pulses  at  a  given  frequency  is 
320  (=1600/5);  hence  the  Nyquist  interval  is  ±160  Hz  and,  for  reasons  to  be  explained,  the  Doppler 
resolution  is  10  Hz.  The  spectra  shown  in  Fig.  5  and  6  were  taken  from  the  backscattered  return  at  2.1 
and  3.0  GHz,  respectively.  The  relation  between  Doppler  frequency  8  and  velocity  v  is  8  =  2v/A.  At 
2.1  GHz,  1  Hz  =  0.1388  knot,  and  at  3.0  GHz,  1  Hz  =  0.0972  knot.  Hence  the  Doppler  resolution  of 
10  Hz  is  approximately  1.4  knots  at  2.1  GHz  and  1.0  knot  at  3.0  GHz. 


The  Doppler  power  spectra  are  computed  according  to 


l  (»)  = 


N 


N- 1 

£  z(nt)  e_2ir“"" 
0 


(2) 


where  N  is  the  number  of  pulses  used  in  the  calculation,  z(t)  =  x(t)  +  iy(t)  is  the  complex  data 
formed  from  the  /-  and  0-channel  voltages  x(t)  and  y(t),  and  t  is  the  interpulse  period  (=3.1  ms). 
Figures  5  and  6  show  smoothed  spectra  for  each  of  five  records  from  the  third  range  cell.  The  smooth¬ 
ing  is  accomplished  by  decomposing  each  record  into  a  number  M  of  smaller  subrecords  containing  N 
pulses  each,  where  A  is  a  power  of  2,  and  averaging  the  A/"periodograms"  (Eq.  (2))  to  obtain  the  spec¬ 
tra  shown  in  the  figures.  In  Figs.  5  and  6,  M  =  17  and  N  =  32,  and  the  Doppler  resolution  is  therefore 
10  Hz.  In  each  case  the  I(v)  are  computed  at  increments  of  v  equal  to  the  resolution.  For  graphical 
display  the  maximum  value  of  I(v)  is  normalized  to  unity. 


When  x(t)  and  y(f)  are  two  stationary  Gaussian  zero-mean  processes  having  the  same  autocorre¬ 
lation  function  and  zero  crosscorrelation,  then,  by  expanding  the  right-hand  side  of  Eq.  (2),  one  can 
easily  establish  that  !{v)  has  the  form  U1  +  F2,  where  U  and  Fare  two  zero-mean  Gaussian  variables 
which  satisfy 

E[U2]  »  £ [ F2] ,  E[UV)  =  0. 

Hence  for  each  value  of  v  the  variable  w  =  l(v)  has  the  exponential  probability  density  function  given 
by  Eq.  (1),  and  in  particular  its  mean  and  standard  deviation  are  equal.  When  the  spectra  are  averaged 
as  described,  the  mean  remains  the  same  but  the  standard  deviation  is  reduced  by  the  factor  1/VA/, 
where  again  M  is  the  number  of  subrecords  into  which  the  original  record  has  been  divided.  The 
averaging  process  therefore  reduces  the  levels  and  numbers  of  spurious  spectra  peaks  but  at  the  cost  of 
reduced  resolution.  These  considerations  suggest  that  some  of  the  spectra  shown  in  Fig.  5  are  truly 
multimodal;  that  is,  at  least  some  of  the  peaks  off  the  main  lobe  are  physically  real  and  not  the  spuri¬ 
ous  effects  of  random  fluctuations.  The  width  of  these  peaks,  when  they  occur,  is  generally  about  20 
Hz,  or  3  knots. 


The  existence  of  multiple  peaks  is  less  apparent  in  the  return  at  3.0  GHz,  and  all  of  the  spikiness 
in  Fig.  6  can  be  reasonably  attributed  to  purely  random  fluctuations.  In  fact  most  of  the  spectra  com¬ 
puted  from  the  3.0-GHz  return  was  even  more  clearly  "unimodal"  than  those  shown  in  Fig.  6.  The 
"most  typical"  spectrum  at  3.0  GHz  is  shown  schematically  in  Fig.  7  (which  was  obtained  by  averaging 
the  spectra  from  the  first  five  records  from  the  sixth  range  cell).  The  hump  appearing  in  the  negative 
frequency  range  is  significant;  for  although  the  radar  was  looking  into  the  wind,  in  every  spectrum  that 
was  examined  (at  either  frequency)  a  significant  portion  of  the  backscattered  signal  energy  lies  in  the 
negative  frequency  range. 

The  wind  shear  can  be  crudely  estimated  by  supposing  the  wind-velocity  null  to  occur  somewhere 
between  the  upper  and  lower  3-dB  points  of  the  radar  beam.  For  the  3.0-GHz  beam  at  the  range  and 
elevation  angle  indicated  for  file  6  in  Table  2,  these  two  points  occur  at  the  heights  1.28  and  2.44  km. 
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Therefore,  given  that  the  surface  wind  speed  was  23  knots  -  12.4  m/s,  the  wind  shear  was  between  5 
and  10  m/s  per  km.  The  extent  of  the  radar  beam  between  the  first  upper  and  lower  nulls  is 
approximately  double  that  between  the  two  3-dB  points,  and,  in  particular,  the  lower  null  of  the  2.1- 
GHz  beam  lay  just  above  the  surface.  It  thus  appears  that  the  vertical  coverage  of  the  beams  contained 
rain  falling  near  the  surface  and  extended  well  into  the  cloud  above. 

The  estimated  value  of  wind  shear  can  also  be  used  to  gauge  the  vertical  extent  A/»  of  the  clutter 
producing  a  spectral  peak  whose  velocity  spread  is  Av:  A  h  —  Av/shear.  For  the  narrow  3-knot  peaks 
discussed,  Ah  —  160  or  320  m  for  the  wind-shear  values  of  5  or  10  m/s  per  km. 

The  five  spectra  shown  in  Fig.  8  were  obtained  from  synthetically  generated  Gaussian  data  with  a 
Gaussian  spectrum  and  are  presented  here  for  comparison  with  the  spectra  shown  in  Figs.  5  and  6. 
Each  of  these  spectra  is  again  an  average  of  17  32-point  periodograms,  and,  as  previously  mentioned,  at 
each  frequency  such  spectra  have  a  standard  deviation  equal  to  the  true  spectral  value  divided  by  VT7. 
The  three  Gaussian  curves  shown  in  each  plot  in  Fig.  8  are:  the  "true"  spectrum,  the  true  spectrum 
plus  one  standard  deviation  unit,  and  the  true  spectrum  minus  one  standard  deviation  unit.  There 
turns  out  to  be  a  15%  probability  that  the  estimated  spectrum  will  exceed  the  true  value  by  at  least  one 
such  unit  and  also,  coincidently,  an  almost  equal  probability  that  the  estimated  spectrum  will  fall  below 
the  true  value  by  at  least  the  same  amount.  It  can  therefore  be  expected  that  70%  of  the  estimated 
spectral  values  will  fall  between  the  upper  and  lower  curves  shown  in  Fig.  8.  The  estimated  spectra  are 
as  spikey  as  expected,  but  their  underlying  Gaussian  character  can  be  discerned. 

Figure  9  shows  two  spectra  taken  from  file  3,  the  single-frequency  file.  Both  spectra  are  computed 
from  the  same  set  of  data,  but  the  first  is  the  average  of  six  64-point  periodograms,  whereas  the  second 
is  the  average  of  13  32-point  periodograms.  The  Doppler  resolutions  are  respectively  25  Hz  (  =  2.5 
knots)  and  50  Hz  (=5  knots),  and  the  Nyquist  interval  has  been  expanded  to  ±800  Hz,  since  the  sin¬ 
gle  frequency  is  transmitted  at  a  PRF  of  1600.  In  all  of  the  spectra  from  this  file,  almost  all  the  signal 
energy  is  in  the  frequency  interval  between  ±200  Hz,  which  increases  our  confidence  that  there  is  no 
appreciable  spectral  foldover  in  the  spectra  shown  in  Fig.  5,  at  a  lower  carrier  frequency.  Due  to 
antenna  rotation,  each  of  the  M  subrecords  into  which  the  total  record  is  divided  originates  from  a 
different  circular  arc  in  space.  (Counting  the  beamwidth,  each  of  these  arcs  is  about  2°  wide,  which  at  a 
range  of  30  km  amounts  to  1  km.)  Therefore,  the  averaging  process  will  tend  to  smooth  out  features 
which  vary  with  azimuth.  The  narrowness  of  the  central  peak  is  a  striking  feature  that  appears  in  many 
of  the  spectra  from  this  file,  again  indicating  that  a  large  fraction  of  the  total  clutter  return  originates 
from  clumps  or  horizontal  bands  whose  vertical  extent  is  rather  small. 


TEMPORAL  AND  SPATIAL  AUTOCORRELATION 


The  power  spectrum  and  autocovariance  of  a  stationary  process  are  a  Fourier-transform  pair. 
However,  the  temporal  autocorrelation  functions  p(m)  plotted  in  Figs.  10  and  11  are  computed  directly 
from  the  data  according  to  the  formula 

p(mt)  -  RJR0, 


where 


Rm  =  - -  T  x\nt]  ■  x[(n  +  m)t]  (3) 

N  -  m  h 

and,  as  before,  x(t)  is  the  /-channel  voltage.  The  autocorrelation  of  the  complex  process 
z(t )  —  x(t)  +  />(/)  should  be  the  same,  on  the  assumption  that  the  /-  and  Q-channel  voltages  are  two 
independent  processes  with  the  same  second-order  statistics. 
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Fig.  6  —  Doppler  power  spectra  at  3.0  GHz.  Each  spectrum 
is  the  average  of  17  32-point  spectra.  The  data  are  from 
records  I  through  S  of  range  cell  3,  file  6. 
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Fig.  9  —  Two  typical  spectra  from  file  3,  the  single-frequency  file,  with  both  spectra  being  calculated  from  the 
same  data  (from  the  same  record  and  range  cell).  The  antenna  was  rotating  at  6  RPM. 


Plots  of  spatial  autocorrelation  functions  are  shown  in  Fig.  12.  These  are  computed  as  by  Eq.  (3), 
but  in  the  present  case  the  lags  are  spatial  rather  than  temporal.  The  compressed  pulse  width  of  0.7  ps 
corresponds  to  a  resolution  cell  width  of  10S.0  m;  however,  as  was  mentioned,  there  is  50%  overlap 
between  successive  range  cells.  Hence  each  spatial  lag  in  these  plots  is  equal  to  52.5  m.  The  plots 
shown  in  Fig.  12  are  smoothed  values  of  p  which  were  obtained  by  averaging  the  results  of  ten 
consecutive  pulses. 

The  most  commonly  accepted  model  for  the  Doppler  power  spectrum  is  the  Gaussian  model,  and 
the  corresponding  correlation  function  [I,  p.  334]  is  a  damped  sinusoid  of  the  form 

pit)  -  g_2(2,r‘r*'/x)  cos  (4jt  V0t/\),  (4) 

whei  y0  is  the  mean  velocity  and  trv  is  the  RMS  velocity  of  the  clutter.  According  to  Nathanson  [1, 
p.  334]  the  value  of  pit)  is  not  very  dependent  on  the  actual  spectrum  shape.  However,  the  values  we 
obtain  for  pit)  are  often  much  smaller  than  the  theoretical  values  calculated  according  to  Eq.  (4)  (with 
<rv  and  V0  computed  from  the  data).  This  has  serious  implications  for  the  predicted  performance  of 
MTI  filters  employing  delay-line  cancelers.  For  example,  the  clutter  attenuation  factors  for  single  and 
double  delay-line  cancelers  are 

CAX  -  [1  -pit)]'1 
and 

CA2  -  (1  -  (4/3)p(f)  +  (I/3)p(2f)]-1, 

and,  without  pressing  very  hard,  we  have  found  that  the  difference  between  the  theoretical  and  actual 
values  of  p  can  easily  produce  a  corresponding  difference  of  5  dB  for  CAX  and  10  dB  for  CA2.  For 
spectra  having  a  variable  number  of  peaks  whose  positions  are  also  variable,  it  might  be  desirable  to  use 
an  adaptive  MTI  to  achieve  clutter  rejection. 

SPATIAL,  TEMPORAL,  AND  FREQUENCY  VARIATION 

Notation 

Three  types  of  averages  appear  in  our  subsequent  discussion,  indicated  by  boldface  type,  single 
bars,  and  double  bars.  To  illustrate,  we  consider  the  three  types  of  "mean  power": 

For  each  of  the  five  frequencies  fh  P,  —  P,(r,  c)  denotes  the  theoretically  "true"  or  ensemble 
mean  power  which  is  characteristic  of  the  return  at  the  ith  frequency  from  the  rth  record  from  the  cth 
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Fig.  12  —  Spatial  autocorrelations.  The  data  are  from  file  1,  and 
each  correlogram  is  a  ten-pulse  average  (pulses  l  through  10). 
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range  cell.  The  sample  mean  P,  —  /*,  ( r,  c)  is  the  numerical  average  of  the  N  values  of  instantaneous 
power  contained  in  the  /th  record  from  the  cth  range  cell,  where  N  —  85  for  files  1  and  2  and  N  *  562 
for  files  5  and  6.  Similarly,  Pt  denotes  the  numerical  average  of  alt  the  values  of  instantaneous  power 
(all  records  from  all  range  cells)  contained  in  a  given  file.  Thus,  each  file  gives  rise  to  Mr  •  Mc  values 
P,(r,  c),  where  1  <  r  <  Mr  and  1  <  c  <  Mc,  but  to  only  one  value  of  Pt: 


P, 


1 

M,  ■  Mc 


11^,  c). 


In  theory  one  attempts  to  relate  the  quantities  P,  M>  the  physical  properties  of  the  clutter.  In  prac¬ 
tice,  the  true  means  P,  are  approximated  by  the  sample  means  Pn  and  the  experimental  verification  of 
any  physical  model  is  complicated  by  random  statistical  fluctuations  in  the  quantities  Pt  which  are  due 
to  the  Gaussian  nature  of  the  voltages  and  would  exist  even  in  the  absence  of  receiver  noise. 

In  general,  data  from  different  files  is  never  collated,  and  the  above  notation  will  be  followed  con¬ 
sistently.  These  general  rules  are_violated  only  in  the  case  of  receiver  noise.  The  quantities  N,  are 
approximated  by  the  file  averages  obtained  from  the  noise  file,  and  we  shall  use  the  same  symbol  N, 
for  both  the  true  value  and  its  estimate  N,.  Thus,  for  example,  we  shall  write  P,  —  N,  instead  of 
P,  -  /V^for  the  "noise  corrected^  value  of  the  indicated  sample  mean.  In  practice,  the  difference 
|N,  -  N,\  is  small,  since  each  /V,  is  the  numerical  average  of  26,400  values  of  instantaneous  noise 
power. 

Nonhomogeneity  of  the  Clutter 

Figure  13  represents  the  spatial  arrangement  of  rain  based  on  the  envelope-detected  radar  return. 
The  data  are  from  file  3,  the  single-frequency  file.  A  standard  lineprinter  has  been  used  to  display  the 
results.  Each  print  line  is  one  radar  pulse,  with  each  character  position  in  the  line  representing  a  range 
cell.  The  data  from  each  range  cell  of  each  pulse  is  quantized  to  a  three-level  scale,  with  the  first  level 
(blank)  consisting  of  returns  below  9,  the  second  level  (plus  sign)  consisting  of  returns  between  9  and 
20,  and  the  third  level  (letter  M)  consisting  of  returns  greater  than  20.  The  data  appear  to  contain  vari¬ 
ous  eddies.  Since  constant-false-alarm-rate  (CFAR)  detectors  assume  a  homogeneous  distribution  in 
the  detector  reference  cells,  rain  clutter  will  produce  a  larger  number  of  false  alarms  than  expected. 
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44  4  4*44  444444  44444***4  ***4  4  444***  4*  44**4***4  4444 

4444  444  4*444  44444****44***4  4  444**1*4  4444  44**44*44  4**44 

4444  4*4  4**44  4  44  4***444**4  44  4*4**44  4444  44**44  44  4*4 

4444  4**44***4  44  44**44444*444*44444444  4*44  4****44*4  *** 

444  4 4 **4***4 4  4  4*44  4  44*4  4**4444  44  *44  4**4**44*4  4**4 

44  4**4**44  444**4  ***** *4  4***4  4*  **  44444**4**4  *** 

444  4*4444  4  44  4******4  4*  ***4  4**444  44444*****44**44 

44  4444444  44  4444***4  4**4  4***4  **444  44444****444**** 

44  444  444  44  444*****4  4**4  ***«4  **444  444  44**4  444*** 

44  44  4444*4  4  4**"***4****4***44  4**444  44  4***44  4  ***4 

44  44  4  4444444**4********4***4  44 ****444  4  4***44444***4 

44  444  444*44  ****************  4***444444  44**44  4***44 

44  4  4  **  4  *K44********"*M*4  4*4444444  4**444  4***4 

44  4  ♦  4*4  4 4** 4**4****** *"**4 44* ♦  4*44  4*4  *  4  ***4 

44  4  4*  44**44**  ******4***44444444*44  4  *444  4*** 

444  4*44**  4**4 *****4  ******4444  444444444*4  44*** 

4  44444  4*44**  4*444****  *******444  444444  44*4  4*  ***4 

4  ******  4 **44**444  4  4***  *4  4**  *4***  4**  4**44  4444444*** 

444  4444*4  44*444  444***4*4*4  ******  4*444  444*4***44 

*44  44  4  ******  4**4444  ****4*444  ****44444444  4*  ***44 

Fig.  13  —  Spalial  power  plot  of  rain 


Table  5  shows  values  of  the  record-averaged  quantities  Pj  -  N3  tabulated  against  range  cell,  for 
each  of  the  two  records  4  and  J  in  files  1  and  2.  A  sharp  gradient  can  be  perceived  in  the  data  from  file 
1.  The  standard  deviation  in_P3  —  N3  produced  by  purely  random  fluctuations  in  these  two  files  is  gen¬ 
erally  of  the  order  of  (0.1)P3,  and  the  variation  in  P3  —  N3  from  record  to  record  and  range  cell  to 
range  cell  is  typically  larger  than  this. 


The  Ideal  Power-Versus-Frequency  Relation 

From  a  straightforward  analysis  of  the  effects  of  receiver  noise  and  calibration  errors,  one  obtains 
the  relation 

P /  -  Z, FW  +  N„  (5) 

where  F  -  F(/)  is  the  "ideal"  power-versus-frequency  relation  that  would  be  observed  in  the  absence  of 
calibration  errors  and  noise  and  Z,  is  the  calibration  error  at  the  Ah  frequency.  Hence,  Z,_=  1  if  the 
calibration  is  perfect,  and  N,  =  0  in  the  absence  of  noise.  Our  problem  is  to  use  the  data  P{  to  obtain 
information  about  the  function  F  -  F(/),  which  depends  on  the  radar  range  R  as  well  as  other  parame¬ 
ters  descriptive  of  the  clutter. 


16 


NRL  REPORT  8639 


Table  5  —  Record- Averaged  Noise-Corrected  Power  ( P3  -  N3) 


Range 

Cell 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


Record  4 


163.93 

158.16 
215.29 
181.95 
168.23 
210.08 
232.01 

259.35 
261.11 
253.45 

200.16 
223.20 
278.63 
296.53 
317.79 
428.91 
455.43 

602.76 

695.36 
682.16 

466.94 
521.89 

803.76 
973.81 

1068.90 

1152.64 
1 100.74 

1079.18 

1017.65 
1112.04 

1489.19 
1769.84 


Record  5 


110.2 

105.21 
121.16 
140.36 

112.73 
108.14 
164.65 
192.03 
246.94 
261.67 
321.06 

319.85 
291.59 
297.09 
232.75 
211.11 
311.56 
378.09 
467.98 

453.22 
540.24 
713.58 

809.74 

778.86 
989.52 
962.83 
796.81 
727.39 
795.35 

1245.36 

1367.94 

1193.32 


Record  4 


73.28 

222.55 
270.30 
254.44 
301.91 
381.97 
384.43 

341.37 
299.08 

353.78 
526.65 
538.09 

523.79 

602.55 
621.46 
672.00 
832.99 

866.55 

639.72 

591.72 
778.02 

831.37 
702.62 

685.38 
792.05 
926.96 
746.35 

706.39 
831.69 


Record  5 


78.76 

199.23 

246.79 
233.01 

278.24 

333.75 
351.03 
346.16 

390.74 
455.12 
430.40 

418.65 
551.30 
586.94 

651.76 
815.69 

840.75 
735.04 

659.79 
861.64 

891.25 
781.01 
799.22 
897.97 

838.75 
655.86 
492.04 
531.99 
498.71 
467.36 

446.66 
460.21 


Record  4 


1620.96 

1488.48 
1334.06 
1631.59 
1646.35 

1445.51 

1273.94 
1173.75 
1091.18 
1135.10 

1135.81 
1128.04 

1289.82 

1289.44 
1732.34 

2118.45 

1660.51 

1354.62 
1531.50 
1765.66 
1693.00 

1545.73 

1411.74 
1386.58 

1384.94 

1389.48 
1465.61 

1552.63 
1447.44 
1420.55 
1507.13 
1611.69 


Record  5 


1052.68 

1140.50 

1526.23 

1629.15 
1275.93 
1189.55 
1122.02 
1302.71 
1554.43 
1538.92 
1378.86 

1393.89 
1583.34 
1632.78 
1445.54 
1314.08 
1263.01 

1188.23 
1234.49 
1438.20 

1626.90 
1459.48 
1478.97 

1617.51 

1645.16 

1402.46 
1300.81 

1389.47 
1237.25 
1167.37 

1072.47 
979.02 


Record  4 


416.80 

391.64 

371.10 
344.59 
367.83 
348.08 

346.94 

443.90 

504.90 
470.01 

341.11 
382.57 

384.74 

325.75 

414.50 

411.51 
348.02 
365.67 
344.20 
297.53 
299.35 
243.01 
296.14 
430.98 
420.86 

481.95 

481.75 
373.71 
349.18 
306.37 
365.45 
427.05 


Record  5 


466.17 

438.13 

362.48 
351.86 
389  25 

369.44 

390.45 
349.97 
383.08 

402.79 
410.01 
353.60 
289.66 
300.09 

347.80 
332.99 
379.02 
361.28 
312.34 
322.77 

358.48 
443.04 

446.76 
400.71 

401.76 
317.12 
283.85 
250.32 
315.91 
324.04 
330.07 
399.82 


Table  6  shows  the  calibrated  values  of  instantaneous  power  averaged  over  all  the  data  in  each  file, 
that  is,  Pf  These  data  are  tabulated  only  for  the  mullifrequency  files,  including  the  noise  file,  file  8. 
The  noise  power  is  not  flat  across  the  frequency  range  because  the  calibration  involved  certain  factors 
(such  as  antenna  gain)  which  do  not  affect  the  receiver  noise.  Also  the  signal-to-noise  ratios  PJ N,  vary 
from  file  to  file.  File  5  is  by  far  the  noisiest  file,  the  signal-to-noise  ratios  in  this  file  all  being  less  than 
9  dB.  The  "relative  values  of  noise-corrected  power"  are  defined  by 

Q,  =  (Pi r-  N,)/(?|  -  N,)  (6a) 

and 

Qi  -  (Pi  -  N,)/(P,  -  N|).  (6b) 

The  quantities  Q,  are  tabulated  in  Table  7.  As  can  be  seen,  the  corresponding  values  of  are  fairly 
constant  from  file  to  file.  This  suggests  the  model 

/»(/,)  -  AQW  +  N„ 
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Table  6  —  File-Averaged  Values  of  Calibrated  Power  P, 


Frequency 

(Ghz) 

Power  P, 

File  1 

File  2 

File  5 

File  6 

File  8 

HHH 

141.5 

77.5 

14.9 

39.7 

451.2 

244.9 

42.5 

127.9 

3.0 

889.5 

480.6 

64.1 

268.5 

12.6 

3.45 

1241.1 

705.1 

95.5 

389.4 

12.9 

3.9 

814.1 

434.8 

79.8 

246.3 

24.3 

Table  7  —  Relative  Values  of  File-Averaged, 
Noise-Corrected  Power 


Frequency 

(GHz) 

* 

Power  Qi 

File  1 

File  2 

File  5 

File  6 

1.00 

1.00 

3.24 

3.26 

3.44 

3.44 

BuHl 

6.49 

6.57 

5.99 

7.66 

3.45 

9.72 

9.61 

11.27 

3.9 

5.84 

5.76 

6.45 

6.64 

where  A  is  a  frequency-independent  parameter  which  varies  from  file  to  file  and  Q  -  £>(/)  is  a  func¬ 
tion  which  depends  on  frequency  only  and  no  other  parameters.  However,  for  reasons  mentioned  in 
the  Introduction,  we  are  not  very  interested  in  the  behavior  of  heavily  averaged  data.  We  therefore 
restrict  our  attention  to  data  contained  in  a  single  file  (or  portion  thereof),  and  we  consider  whether  the 
data  within  a  given  file  conforms  to  the  model 

P,--4Q, +  N„  (7) 

where  the  parameter  A  varies  with  time  and  position  but  not  with  frequency  whereas  the  function 

Q,-Q(y;)-(p,  -n,)/(p,-n1)  (8) 

depends  only  on  frequency.  Comparing  Eq.  (7)  to  Eq.  (5),  we  see  that  the  model  given  by  Eq.  (7) 
asserts  that  the  "ideal"  mean-power- versus-frequency  relation  is 

F(/)  -  A  •  (Q ,/Z,).  (9) 

where  the  first  factor  contains  all  the  spatial  and  temporal  dependence  and  the  second  factor  contains  all 
the  frequency  dependence. 

Parenthetically,  the  validity  of  the  model  given  by  Eq.  (7)  was  somewhat  unexpected.  In  the 
literature  one  finds  discussions  of  various  scattering  mechanisms  (discussed  in  the  next  subsection), 
each  of  which  give  rise  to  an  ideal  law  of  the  type  P  -  A,r  (in  the  absence  of  noise  and  calibration 
errors),  with  different  effects  giving  rise  to  different  values  of  n.  If,  say,  two  such  effects  are  operating 
simultaneously,  then 

p,  -  Aur+  A2r  +  n„ 

which  does  not  conform  to  Eq.  (7),  since  the  noise-independent  part  cannot  be  expressed  as  the  pro¬ 
duct  A  •  Q(f),  where  A  is  frequency  independent.  Before  examining  the  data,  we  expected  to  observe 
a  law  of  this  type  with  m  -  4  and  0  ^  n  <  2.  Physically  this  corresponds  to  certain  turbulence  effects 
superimposed  upon  a  more-or-less  uniform  distribution  of  clutter  in  the  beam. 
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Again,  P,  =  P,(r,  c)  varies  with  record  rand  range  cell  c,  and  the  hypothesis  now  to  be  tested  is 
that,  for  data  wkhin  a  given  file,  Q,  does  not  vary  with  r  and  c.  The  quantities  Q^are  approximated  by 
the  quantities  Qh  and  to  test  the  hypothesis,  the  standard  deviations  of  the  Q ,  are  first  calculated 
theoretically,  on  the  assumption  that  the  Q,  are  constant.  Then  these  theoretically  calculated  standard 
deviations  are  compared  with  the  actual  sample  standard  deviations  computed  directly  from  the  A/,  •  Mc 
samples  of  Q ,  contained  in  the  given  file.  The  results  are  shown  in  Table  8.  In  our  theoretical  calcula¬ 
tions  we  neglected  the  effects  of  receiver  noise,  which  explains  why  the  agreement  is  not  as  good  in  the 
noisy  files  5  and  6  as  in  files  1  and  2.  Keeping  in  mind  the  nature  and  effects  of  these  simplifying 
assumptions,  we  judge  that  the  results  presented  in  Table  8  represents  good  fit  between  the  data  and 
the  model.  Details  of  the  statistical  analyses  are  given  in  Ref.  4.  The  validity  of  this  result  is  not 
affected  by  the  presence  of  the  calibration  errors  Z,,  since  they  are  constants  and  their  effects  are  multi¬ 
plicative.  In  fact  it  can  be  shown  (4]  that  ratios  of  the  calculated  to  the  sample  standard  deviations 
shown  in  Table  8  are  mathematically  identical  to  the  values  that  are  obtained  by  using  the  raw  uncali¬ 
brated  data. 


Table  8  —  Comparison  of  Theoretically  Calculated  Standard  Deviations  of  Q, 
with  Sample  Standard  Deviations  of 


Frequency 

(GHz) 

Standard  Deviation 

File  1 

File  2 

File  5 

File  6 

Sample 

Theo¬ 

retical 

Ratio 

Sample 

Theo¬ 

retical 

Ratio 

Sample 

Theo¬ 

retical 

Ratio 

Sample 

Theo¬ 

retical 

Ratio 

mm 

1.00 

1.86 

2.49 

1  1.70 

0.90 

1.82 

2.52 

1.64 

m 

0.86 

1.94 

2.71 

1.57 

0.87 

1.92 

2.78 

1.67 

0.99 

1-01 

0.98 

0.94 

0.75 

0.96 

1.65 

1.02 

1.70 

1.22 

1.28 

1.23 

■ 

Evidence  for  the  Vertical-Stratification  Model 

As  mentioned  in  the  Introduction,  two  physical  mechanisms  could  explain  the  existence  of  multi¬ 
ple  peaks  in  the  clutter  Doppler  spectra:  wind  shear  operating  on  a  vertically  stratified  medium,  and 
Bragg-type  scattering  from  turbulent  eddies.  We  will  now  briefly  describe  why  we  believe  the  data  favor 
the  vertical  stratification  model.  Details  of  the  analysis,  which  is  rather  prolix,  are  given  in  Ref.  5. 

We  have  shown  that  the  data  within  any  given  file  is  consistent  with  the  operation  of  a  simple 
power- versus-frequency  relation  (Eq.  (9))  involving  only  a  single  multiplicative,  frequency- 
independent,  and  variable  parameter  A.  For  a  perfectly  calibrated  system  the  power- versus-frequency 
law  given  by  Eq.  (9)  assumes  the  form 

F(/)  —  A  •  Q(/),  (10) 

where,  again,  the  function  Q  contains  all  of  the  frequency  dependence.  The  simplest  law  of  this  kind  is 
given  by 

f  (f)-A-.r,  (id 

where  the  exponent  n  is  fixed.  In  the  literature  one  finds  discussions  of  various  scattering  mechanisms, 
each  of  which  give  rise  to  an  ideal  power-versus-frequency  law  of  this  type.  For  the  classical  case  of 
incoherent  scattering  from  drops  uniformly  distributed  throughout  the  radar  resolution  cell,  we  have 
n  —  4  [5].  For  incoherent  scattering  from  narrow  horizontal  bands  or  clumps  of  clutter,  it  turns  out 
that  n  —  5  [4],  whereas  the  Bragg-scattering  effects  discussed  in  the  literature  are  expected  to  produce 
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values  of  n  which  lie  in  the  range  0  <  n  <  2  [6-8],  (In  addition,  certain  coherent  scattering  effects 
discussed  by  Kerr  [3]  would  also  produce,  exponents  in  the  range  0  <  n  <  2.)  In  all  these  theoretical 
formulations,  it  is  assumed  that  the  radar  cross  section  c r  of  each  drop  follows  the  Rayleigh  scattering 
law,  <r  a  /*,  and  due  account  is  taken  of  the  variation  of  the  size  of  the  radar  resolution  cell  with  fre¬ 
quency. 


The  object  of  our  analysis  is  to  decide  between  the  two  alternative  models:  Bragg  scattering  (0  < 
n  <  2),  and  vertical  stratification  of  the  clutter  (n  =  5).  From  the  data  in  each  record  from  each  range 
cell,  we  formed  four  estimates  of  n  according  to 

«a/«  I*/*"  log  (Pi  ~  N,),  (12) 

<-] 


for  M  =  2,  3,  4,  and  5;  that  is,  we  formed  each  estimate  hM  by  using  only  the  first  M  frequencies.  We 
chose  the  weights  so  that  Eq.  (12)  gives  the  correct  result  (hM  -  n)  in  the  absence  of  calibration 
errors  and  receiver  noise,  and  among  all  such  estimates  we  chose  the  one  which  minimizes  the  effects 
of  calibration  errors  (since  the  effects  of  receiver  noise  are  smaller  and  can  be  averaged  out).  It  turns 
out  that  these  "optimal"  weights  are  given  by 


f = 


loe  -  -h  | 108  f' 


^  (log  fj)2  - 

j- 1 


M 


I  M 

Z  !°g  fj 

j- 1 


(13) 


Each  data  file  gives  rise  to  Mr  •  Mc  values  of  hM  (one  from  each  record  from  each  range  cell). 
The  file-averaged  quantities  nM  are  shown  in  Table  9.  For  each  fixed  M  the  file-to-file  variation  of  nM 
is  rather  small,  but  within  each  given  file  the  variation  nM  with  M  is  rather  large  and  systematic.  The 
sample  standard  deviations  in  hM  are  shown  in  Table  10.  The  values  of  these  quantities  fall  within  the 
range  of  statistical  expectation  on  the  assumption  that  the  ideal  power-versus-frequency  law  is  given  by 
F(f)  —  Af",  where  A  is  frequency  independent. 

Table  9  —  Values  of  Exponent  Estimates  nM 


File 

nl 

n) 

n 4 

ni 

1 

6.0 

5.2 

4.5 

3.1 

2 

6.2 

5.4 

4.6 

3.2 

5 

5.8 

4.6 

4.2 

3.0 

6 

6.3 

5.7 

4.9 

3.4 

•i 


Table  10  —  Sample  Standard  Deviations  of  hM 


File 

SD(n2 ) 

SD(hi) 

SD(n4) 

SD(n5) 

l 

1.33 

0.74 

0.49 

0.37 

2 

1.15 

0.68 

0.47 

0.34 

0.79 

0.45 

0.41 

flP-1 

0.38 

0.25 

0.18 

We  now  cowe  to  the  heart  of  our  analysis.  It  can  be  shown  that  the  effects  of  receiver  noise  on 
nM  are  small  and  that  the  errors  nM  -  n  are  almost  entirely  due  to  the  calibration  errors  Zt: 

»*-»  +  £  log  Z, 

I- 1 
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The  second  term  on  the  right-hand  side  has  the  form  of  a  vector  dot  product  between  the  vectors 
and  (log  Z,},  and  using  a  certain  generalization  of  the  Schwarz  inequality,  one  can  obtain  lower 
bounds  to  the  RMS  calibration  errors  (dB), 


RMS  = 


17 1  no  iogio  z,): 


1/2 


(14) 


in  terms  of  the  quantities  nM  —  n.  These  lower  bounds  are  shown  in  Table  11,  with  the  data  being 
from  file  6.  In  this  table  the  assumed  "true”  value  of  n  is  increased  in  integral  steps  from  n  —  0  to 
n  =  10,  and  each  entry  shows  a  lower  bound  to  the  RMS  calibration  errors  (dB)  defined  by  Eq.  (14). 
An  RMS  value  above  2  dB  is  considered  to  be  beyond  the  range  of  credibility,  and  it  is  on  this  basis 
that  we  reject  the  Bragg-scattering  model,  since  values  of  n  in  the  range  0  <  n  ^  2  would  imply  RMS 
calibration  errors  of  3  dB  and  more.  Also,  the  smallest  lower  bounds  are  obtained  near  n  =  5,  which 
corresponds  to  vertical  stratification  (scattering  from  narrow  horizontal  bands  of  dense  clutter).  Finally, 
the  evidence  is  especially  compelling  if  we  exclude  the  line  at  M  =  5,  which  contains  the  results  involv¬ 
ing  the  anomaly  at  the  highest  frequency  (3.9  GHz). 


Table  1 1  —  Lower  Bounds  to  the  RMS  Calibration  Error 


~  £(10  1og10Z,); 


for  Assumed  Exponents  n.  The  data  are  from  file  6. 


1/2 


M 

Lower  Bound  (dB) 

n=0 

n=  1 

n  =  2 

n  =  3 

n=4 

n  =  5 

n  =  6 

n  =  7 

n  =  8 

n“9 

n=  10 

B9 

2.7 

2.2 

1.8 

1.4 

mm 

K29 

0.1 

0.3 

1.1 

1.6 

u 

4.4 

3.7 

2.9 

mm 

0.4 

1.1 

KX 

3.4 

n 

5.6 

4.5 

3.4 

■IS 

1.4 

5.8 

H 

5.3 

3.8 

2.5 

2.0 

2.1 

4.1 

7.0 

10.0 

SUMMARY  OF  THE  MAIN  CONCLUSIONS 

The  conclusions  to  be  listed  are  not  intended  to  represent  general  facts  which  hold  in  every  case, 
and  one  would  expect  that  measurements  taken  under  different  meteorological  conditions  could  yield 
different  results.  Moreover,  besides  these  conclusions  themselves  it  is  of  equal  importance  to 
emphasize  the  nature  of  the  statistical  methods  used  to  obtain  them,  which  were  designed  to  be  valid  in 
the  presence  of  calibration  errors  and  without  the  assumption  that  the  clutter  process  was  stationary 
(homogeneous)  over  extended  regions  of  space  or  durations  of  time.  In  particular,  we  did  not  attempt 
to  establish  trends  or  functional  relationships  by  "fitting  curves"  to  the  data.  Even  in  the  absence  of 
calibration  errors,  such  regression  techniques  could  establish  only  the  "average"  behavior  of  the  clutter, 
which  would  be  of  little  use  either  in  gaining  an  understanding  of  the  rapidly  changing  clutter  process  or 
in  designing  radar  systems  which  must  work  in  a  highly  variable  and  complex  environment.  For  exam¬ 
ple,  the  heavily  averaged  spectrum  shown  in  Fig.  7  indicates  the  presence  of  a  wind  shear  but  gives  no 
hint  to  the  vertical  stratification  suggested  by  the  spikiness  of  some  of  the  spectra  shown  in  Fig.  5. 

Of  the  four  points  of  summary  that  follow,  the  first  two  (which  are  purely  empirical)  have  poten¬ 
tially  serious  implications  for  the  performance  of  CFAR  and  MTI  detectors,  and  the  last  two  are  con¬ 
cerned  with  determining  what  physical  mechanisms  might  be  involved  in  the  production  of  multipeaked 
clutter  spectra: 

•  The  clutter  was  nonhomogeneous  in  space  and  nonstationary  in  time.  Large  qualitative 
differences  could  sometimes  be  discerned  in  clutter  spectra  which  were  as  little  as  2  s  apart  (which  was 
the  smallest  elapsed  time  between  consecutive  data  records). 
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•  The  clutter  Doppler  spectra  were  non-Gaussian  in  shape,  and  the  spectra  were  frequently  mul- 
tipeaked.  Although  the  radar  was  looking  into  the  wind  at  the  surface,  every  single  spectrum  showed  a 
significant  amount  of  signal  energy  in  the  negative  Doppler  range,  indicating  the  presence  of  a  wind 
shear. 


•  An  analysis  of  .ne  variation  of  backscattered  power  with  position,  time,  and  frequency  revealed 
that  the  data  were  consistent  with  the  operation  of  a  simple  ideal  mean-power- versus-frequency  relation 
of  the  type  P  =  A  ■  Q(f),  where  A  is  a  frequency-independent  parameter  which  varies  with  time  and 
position  and  Q  =  Q(f)  is  some  function  which  depends  only  on  frequency.  (This  result  is  valid  in  the 
presence  of  calibration  errors  and  could  even  have  been  obtained  from  the  raw  uncalibrated  data.) 

•  The  multiple  peaks  which  were  observed  in  the  clutter  spectra  could  have  been  caused  by  the 
wind  shear  operating  on  a  vertically  stratified  medium  or  by  Bragg  scattering  from  turbulent  eddies.  For 
Bragg  scattering,  Q(f)  =  /",  where  n  falls  in  the  range  0  ^  n  <  2,  whereas  scattering  from  horizontal 
bands  or  clumps  of  clutter  would  produce  the  exponential  value  n  =  5.  To  establish  a  reasonable  range 
of  exponential  values,  lower  bounds  to  the  RMS  calibration  errors  were  calculated  in  terms  of  the  quan- 

f  tities  h  —  n,  where  n  is  an  assumed  "true"  value  (varied  from  0  to  10  in  integral  steps)  and  Pi  is  its  esti¬ 

mate  (formed  from  certain  weighted  logarithmic  sums).  Values  of  n  which  implied  RMS  calibration 
errors  larger  than  2  dB  were  rejected  as  reasonable  possibilities,  and  on  this  basis  the  vertical- 
stratification  model  (n  =  5)  is  favored  over  the  Bragg-scattering  model  (0  <  n  <  2),  since  the  latter 
implies  RMS  calibration  errors  larger  than  3  dB.  In  fact,  the  "reasonable"  range  of  exponential  values 
was  found  to  lie  in  the  interval  3‘  <  n  <  7,  and  the  smallest  lower  bounds  to  the  RMS  calibration  error 
occurred  at  values  of  n  close  to  5. 
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